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Abstract 

We calculate the spectrum of Lyapunov exponents for a point particle mov- 
ing in a random array of fixed hard disk or hard sphere scatterers, i.e. the 
disordered Lorentz gas, in a generic nonequilibrium situation. In a large sys- 
tem which is finite in at least some directions, and with absorbing boundary 
conditions, the moving particle escapes the system with probability one. How- 
ever, there is a set of zero Lebesgue measure of initial phase points for the 
moving particle, such that escape never occurs. Typically, this set of points 
forms a fractal repeller, and the Lyapunov spectrum is calculated here for 
trajectories on this repeller. For this calculation, we need the solution of the 
recently introduced extended Boltzmann equation for the nonequilibrium dis- 



1 



tribution of the radius of curvature matrix and the solution of the standard 

Boltzmann equation. The escape-rate formahsm then gives an exphcit result 

for the Kolmogorov Sinai entropy on the repeller. 
PACS numbers: 05.20.Dd, 05.40.-a, 05.45, 05.60 
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I. INTRODUCTION 



In this paper we extend the analysis of the chaotic properties of dilute, random Lorentz 
gases given in (denoted by (I)) to include open systems with absorbing boundaries. The 
Lorentz gas consists of a point particle moving in a system of identical hard disk (d=2) 
or hard sphere (d=3) scatterers of radius a. In a dilute, random Lorentz gas, the average 
distance between the scatterers is large compared to their radius, a, and the scatterers are 
placed at random on the plane or in space without overlapping each other. The interest 
in open systems, with absorbing boundaries, is occasioned by the escape-rate method of 
Gaspard and Nicolis which relates the coefficient of diffusion for the moving particle in 
the Lorentz gas to the dynamical properties of particles on the set of trajectories that never 
escape from the system. 

The method of Gaspard and Nicolis is based on an identity in the theory of open, 
hyperbolic dynamical systems, called the escape-rate formula. This formula is an expression 
for the rate of decay of the probability, -P(t), of finding a moving particle in an open region, 
V ^ surrounded by absorbing boundaries, at time t. If the motion of the moving particle is 
hyperbolic, the time dependence of P{t) is exponential, decaying as exp(— 7/!:), where the 
escape-rate, 7, is given by 

7=5: \{n)~hKsin). (1) 

Ai>0 

Here IZ denotes the set of initial phase points for the moving particle which are on trajectories 
that never escape from the system. This set of points is called a "repeller" in the phase space, 
typically of measure zero with respect to the usual Lebesgue measure, with unstable and 
stable manifolds characterized by Lyapunov exponents A, (7?.), and Kolmogorov-Sinai (KS) 
entropy hxsiJZ)- The sum in Eq. (|l|) is only over the positive Lyapunov exponents, and the 
Lyapunov exponents and KS entropy are to be calculated with respect to an appropriate 
measure on the repeller . This equation may be considered as the microscopic expression 
for the escape-rate of the particle from the open region, V. A macroscopic expression for the 
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escape-rate is provided by the diffusion equation satisfied by P{t) on large time and large 
space scales, 

^ = DV'P{t), (2) 

where D is the macroscopic diffusion coefficient for the moving particle in V. The solution 
of this equation, for long times, and for absorbing boundary conditions is of the form 

P(t) = exp[-i^)Dt], (3) 

where L is a length characterizing the distance to the absorbing boundary of interior points 
in V, and c is a numerical factor determined by the shape of V and the absorbing boundary 
conditions. Since the microscopic and the macroscopic expressions for the escape-rate de- 
scribe the same escape process, they have to be identical, which leads to the Gaspard-Nicolis 
formula 



D = lim 

L^oo c 



J2 Hn)-hKs{n) 

Ai>0 



(4) 



Here the limit L — > oo is taken to remove finite size corrections, including the effects of 
microscopic boundary layers. It is worth mentioning that the Green-Kubo microscopic ex- 
pressions for the diffusion coefficient, and Sinai's expression for the sum of positive 
Lyapunov exponents represent these quantities as infinite time integrals over appropriate 
dynamical quantities. In order to apply Eqs. without having to worry about sub- 

tleties due to the possible slow convergence of these integrals (due to long-time-tail effects) 
we assume that all of the dynamical quantities, Xi,hKs reach their long-time asymptotic 
values on a time scale which is short compared to the hydrodynamic time scale on which 
the diffusion equation applieJil. The escape-rate formula for diffusion, Eq. (|), and the 
generalizations to other transport coefficients show a striking connection between the 
macroscopic quantities that control hydrodynamic processes, the transport coefficients, and 



^This assumption is certainly a reasonable one for tlie low density cases we consider here. 
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the microscopic quantities that describe the chaotic dynamics taking place on the repeller. 
A detailed discussion of this can be found elsewhere 

The purpose of this paper is to provide an analytical calculation of the positive Lyapunov 
exponents on the repeller, Xi(7V), for the random, dilute Lorentz gas, and to use these, 
together with known values for the diffusion coefficient, D, to determine h^siT^)) the KS 
entropy of the trajectories on the repeller. We take the same approach as in I where we used 
kinetic theory arguments to calculate the Lyapunov spectra of two and three dimensional 
random, dilute Lorentz gases in equilibrium. Here we are concerned with a nonequilibrium 
situation where particles escape from the system. We will see that spatial inhomogeneities 
introduced by the absorbing boundaries will require some significant modifications to our 
previous calculations. There we used mean-free-path arguments and some results from the 
theory of products of random matrices to determine low density values for the individual 
Lyapunov exponents and we used an extended Lorentz-Boltzmann equation as an efficient 
method for determining the sum of the positive Lyapunov exponents. Here we will do the 
same for the Lyapunov exponents on the repeller. We mention that before we developed this 
method |Q] analytical results for Lyapunov exponents on repellers had only been obtained 
for simple one dimensional models [^,0 • Otherwise one had to use numerical methods |]TT] . 



It will be helpful to recall some ideas from I. There we obtained the individual Lyapunov 
exponents as well as their sums in terms of various averages over functions of a radius of 
curvature (ROC) matrix p, using an appropriate distribution function. A central notion 
introduced in I is the use of an extended Lorentz-Boltzmann equation (ELBE) to determine 
the distribution of the elements of the radius of curvature matrices needed for our calcula- 
tions. The ELBE was derived heuristically in I. In the absence of external fields acting on 
the moving particle it is given by 



dF ^ „^ . dF 



j dn\v ■ n\ Q{v ■ h) J dp'5{p — p'{p))F{f, vf, p', t) — Q{—v ■ n)F{f, v, p, t) 



■ (5) 



The notation is exactly the same as in I. The primed variables denote the restituting values. 



which lead to the unprimed values after a collision, and n denotes a unit vector from the 
center of a scatterer to the point of impact of the moving particle at a collision. The solution 
of Eq. d^) is normalized according to 

Jdp'F{f,v,p',t) = fB{f,v,t), (6) 



where fs is the solution of the standard Lorentz-Boltzmann equation |T^. In this paper 
we will assume, as in I, that the elements of the inverse of the restituting matrix, 
are typically small compared to the inverse of the scatterer radius a~^. As a result, we may 
simplify the delta function appearing in the collision integral on the right hand side of Eq. 

(!)■ 

The plan of the paper is as follows: In Section II we discuss the problem of calculating 
averages on the repeller in phase space. There we argue that it is necessary to introduce 
a survival probability for a particle that is inside the system at time t, to be still inside 
the system at time t + T. We argue that this probability is needed to guarantee that the 
properties we calculate are actually those of the repeller and not merely some asymptotic 
properties of a set of particles that eventually escape from the region A. (Not taking this 
survival probability into account led to an erroneous result in a previous publication on this 
subject, subsequently corrected in an erratum [§].) We will set up the formalism that will 
allow us to calculate the Lyapunov spectrum on the repeller in Section III. In Section IV, 
we will treat the two dimensional case, and in Section V, the three dimensional case. We 
conclude in Section VI with a summary of our results, a number of remarks, and a discussion 
of some interesting open questions. 

II. AVERAGING ON THE REPELLER 

To treat the open system correctly it is necessary to develop a tool which guarantees 
that the quantities we calculate are the actual Lyapunov exponents for trajectories on the 
repeller. In general, the measure on a repeller is a very singular object. For large systems 
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we expect that this measure will be very similar for all typical configurations of scatterers, 
when observed on length scales small compared to typical macroscopic length scales (system 
size) but large compared to the mean free path between collisions. Therefore averages on the 
repeller may be replaced by averages over smooth reduced distribution functions for the light 
particle alone obtained by averaging over all configurations of scatterers. Further support of 
this picture is provided by the observation that for large systems the fractal dimension is close 



to the embedding dimension, as follows from Young's formula [1^ relating the Lyapunov 
exponents to the information dimension of the repeller. 

For dilute Lorentz gases without escape, the probability density of the ROC matrix 
elements is obtained as the time independent solution of the generalized Lorentz-Boltzmann 
equation for F{p) given by Eq. (^). For a system with escape, the solution of the time 
dependent Lorentz-Boltzmann equation only determines the probability of finding a particle 
at point r with velocity v with a ROC matrix p at time t. It does not contain any information 
about the future behavior of this particle. More specifically it does not exclude the possibility 
of this particle leaving the system at a later time t + r. To obtain the smoothed density on 
the repeller we therefore have to weight this Lorentz-Boltzmann density with the survival 
probability. This is the conditional probability S{r,v,t\t + T), that a particle at point 
r with velocity v at time t will still be in the system at the time t + T . It does not 
depend on p because all trajectories in a bundle are infinitesimally close to each other. 
The introduction of the survival probability in order to obtain a proper description of the 
repeller is very reminiscent of, and essentially identical to, methods used to compute fractal 
dimensions and other properties of attractors and repellers in more traditional dynamical 
systems calculations [p!0| ,p!^. 

The survival probability that we need, ^(r, t|t + T), can, in d dimensions, be written 
as an integral over the conditional probability S{r, v, t\f/, v/, t + T) as 

S{f, V, t\t + T)= J d'^ffd'^vfS{f, V, t\ff, vf, t + T), (7) 

with initial condition 



S{f,v,t\r ,tf ,t) = S{r — r')S{v — v) (8) 

The average of a function g{p) on the repeller is then given by 

/ dpdfdvF{f, V, p, t)S{f, V, t\t + T)g{p) 
{9{p))Rep = /hn ^irn j — — — — (9) 



The hmit T — * oo has to be taken first, to guarantee that only the trajectories which 
never leave the system are counted. The limits of numerator and denominator vanish sep- 
arately, since the repeller is a set of Lebesgue measure zero. The conditional probability 
S{f, V, t\r', v' ,t + T) is the solution of the ordinary Lorentz-Boltzmann equation for open sys- 
tems fs{^) v', t + T) with the initial condition specified in Eq. (j^). Due to time reversibility 
and time translation invariance, this probability is the same as that for finding a particle 
at point f with velocity —v at time t, that was at point P with velocity —v' at time t — T. 
Thus, 

S{f, V, t\7^, v',t + T) = S{^, -tT, t - T\f, -V, t) = S{f, -{T, 0|f, -v, T). (10) 

The integration in Eq. with respect to and allows us to replace S{r, v, t\t + T) by 
the solution, fs{r, —v,T)Af6{v — Vo), of the Lorentz-Boltzmann equation with homogeneous 
initial condition fs{r, —v, 0) = 1. Here A/" is a normalization constant given by A/" = (27rfo)"^ 
and {AttVq)^^ in two and three dimensions, respectively. That is, we need the long time 
solution of the Lorentz-Boltzmann equation, with absorbing boundary conditions, and with 
an initial condition which is uniform in space and in velocity directions. Such a solution will 
not stay uniform due to the escape of particles through the absorbing boundaries. 

Therefore Eq. is equivalent to 

/ dpdrdvF{r, v, p, t)fs{r, -v, T)g{p) 

{dip)) Rep = lim Jim ^ (11) 

/ dpdrdvF{f,v, p,t)fsir,-v,T) 



Eq. (|TT]) reduces the calculation of the sum of the Lyapunov exponents and the maxi- 
mum Lyapunov exponent on the repeller to standard integrations, when the solution of the 
Lorentz-Boltzmann equation for fs{r, —v, T) and that of the ELBE for F are obtained. In 



the following sections we will solve these equations for large systems with absorbing bound- 
ary conditions. For long times, t, the solution of the ELBE for F will be obtained as a 
generalized Chapman-Enskog [|T^ expansion of the form 
F{f, V, p, t) = J\f5{v - VQ){ipQ{p)n^{r, t) + V^i(p){/ ■ Vn„(f, t) + V^2(p)^^V2n^(r, t) + ■ ■ ■)• 

(12) 

Also nm{r,t) is the slowest decaying eigenmode of the diffusion equation with the given 
absorbing boundary conditions, 

UmidV) = 0. (13) 

In second order in the gradient we have kept only the scalar part V^n^. A possible contri- 
bution of order to the Lyapunov exponents from a term in Eq. ([T^ ) proportional to the 
traceless tensor vv — {v'^/d) 1 vanishes after integration with respect to the velocity and is 
therefore neglected in this equation. The solution of the usual Lorentz-Boltzmann equation 
for fsir, —v,T) can also be written as a Chapman-Enskog expansion 

fs{f, -V, T) = M5{v - vo){n^{f, T) - CdV ■ Vn™(f, T) + ■■■). (14) 



The constants in Eq. ([T^ ) are C2 = — 3/(4z/) and C3 = —l/i> in two and three dimensions 
respectively, where v is the mean collision frequency. Due to the normalization condition, 
Eq. (^, the functions ipi in Eq. (|T2|) have to fulfill the conditions 



dpVo(p') = 1, (15) 
Jdp'Mp')=Cd, (16) 
I dp'Mp') = 0. (17) 

It is important to note that in Eqs. (|T2|) and (|14|), the eigenmodes nm{r, t) and nm{r, T) 
have the functional forms uoij) exp[— tcu] and no(r) exp[— Tcu], respectively, where uj denotes 
the eigenvalue of the slowest decaying mode of the diffusion equation, and no(r) is the 
corresponding eigenfunction. 
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Using Eqs. (|T2| - we can write Eq. (0) as 

[dp [ d'r kn^lr-) - vH^^P^ + V^2)(Vno(r-))^ 

(gip)) ^ L d 



9{P) 

+ ■••. (18) 



Rep 2 2 

j d'T{nl{f)-'^{Vno{f)f) 

We note that exponentially decaying factors have canceled in the numerator and denominator 
ofEq. (|T|). 

Since we only kept terms up to the order when deriving (plSf ), we also have to expand 
the denominator in Eq. (|18|) . The final result for averaging a quantity g{p) on the repeller, 
up to and including terms of order V^, is given by 



{dip)) Rep = J dpyjQ + 

with 



d Cd 



v'of} 9{P) + ---, (19) 



c/V(Vno(f))' 

t = ' f , , • (20) 

/ d'^rnlif) 

The quantity q may be interpreted as a wave vector. If, for example, a system is considered 



with absorbing boundaries at x = ±L/2, but of infinite extent in other directions, Eq. (20) 
can be evaluated to show that q = g™" = {t^ / L) is the smallest possible wave vector in the 
X direction. We will consistently neglect effects due to microscopic, kinetic boundary layers 
near dV , since such effects are unimportant in the escape-rate formalism as the size of the 
system becomes large. 

That this is the case can be understood in the following way: The fraction of the volume 
taken up by the boundary layers is of order \/L with A the mean free path and L a length 
of the order of the diameter of the system. The density in the boundary layer is equally 
of order \/L compared to the average density in the system, as a result of the absorbing 
boundary condition. Finally the escape rate near the boundary is of similar order, as it is 
proportional to the local density. As a result the effects of the boundary layer are of order 
(A/L)^, whereas here we will be interested in terms of order (A/L)^ only. 



10 



III. THE LYAPUNOV EXPONENTS 



The strategy presented so far is applicable to the calculation of all quantities which can 
be written as ensemble averages of functions of the ROC matrix. In our case this always 
happens, if a quantity is given as a time average over a function of p{t). The sum of the 
positive Lyapunov exponents fulfills this property as shown by Sinai |T^. We assume that 
trajectories on the repeller for the Lorentz gas are sufficiently ergodic, so that we can write 
the sum of the exponents as an average over the appropriate ensemble of ROC matrices (see 
also I). 

A,(7^) = v{Trace{p-'))R,p. (21) 

Ai>0 

In two dimensions this is trivially also the largest exponent. The maximum Lyapunov 
exponent in 3 dimensions is not calculated as a time average, but as an average over a 
function of time of free flight and collision parameters. The separation of trajectories at a 
given time can then be written as product of matrices, each describing the propagation of 
the ROC between two collisions, as explained in I [Q. Using an identical method, we find 

N 

Sf^(t) = n Ui{Ti, (Pi, ai)5r^{0), with 

i=l 

Ui{T,, 0„ a,) = Texp( vp-\t', 0„ a,)dt'). (22) 

Jti 

Here (pi and are the scattering angles at the i-th collision and Xj = tj+i — is the time 
of free flight between the i-th and i + 1-th collision. It is important to notice that Eq. (|22D 
cannot be written in terms of a time integral over a local function of p{t), since the Ui 
matrices typically do not commute with each other. 

Since for dilute Lorentz gases correlations between collision events are not important 
in the limit of low densities, we can make the approximation that the matrices Ui are 
independent, randomly, and to the order in gradients that we use here, isotropically dis- 
tributed matrices, each with independent free flight times and collision parameters selected 



from appropriate distributions, to be discussed below [|15[. We will postpone a discussion 
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of the isotropy of the U matrices until Section V and Appendix B, where we will calculate 
the largest Lyapunov exponent for three dimensional systems. The maximum Lyapunov 
exponent is then given by 

1 ^ 

A^a.(7^) = Ihn ^^^ln(|t/,-el), (23) 

where e is an arbitrary unit vector normal to the velocity v. As in the case of infinite 
systems, we now write the right hand side of Eq. (^) as an average over the distribution 
of the matrices L/j. To do this we have first to derive the distribution of particles, which 
collided a time of free flight r ago, v, r, t), for an open system. This is done in Appendix 
A. Then we have to make sure that the average is restricted to the repeller. This can be 
achieved as described in Section II by using the survival probability. We argue as follows: 
To determine the appropriate average of ln(|Uj ■ e|), one needs the distribution of particles 
that have collided at a point r, with the scattering angles 0, a, that have the velocity v 
after collision, and that travel freely for a time r until the next collision. Further, to include 
only trajectories on the repeller, we again need the survival probability, ^(r, v, t\t + T), for a 
particle with r, v at time t to remain in the system at least until time t + T. We can express 
the distribution of particles colliding at point r and having velocity v after collision in terms 
of the distribution of particles arriving to collide at point f with the restituting velocity v/, 
since the rate at which particles arrive at a scatterer before collision should be equal to the 
rate at which they leave the scatterer after collision. The former rate can be calculated from 
kinetic theory without difficulty due to the assumption of molecular chaos, while the latter 
rate requires that correlations between particles and scatterers, produced by a collision, be 
taken into account. Note also that I'll ■ n| = \v/ ■ n\, and that v ■ "h > 0, while v/ ■ n < 0. 
We finally arrive at the following general expression for the maximum Lyapunov exponent 
of the open system, 

Amaa;(^) = '^cr(ln \U ■ e\) cr 

/ drdhdfdvln ■ v\Q{n ■ v)fF{r, vf, r, t)fs{f, —v, T)\n\U ■ e\ 

lim Ucr- J . (24) 

dTdndfdv\n ■ v\Q{n ■ v)fp{f, vf, r, t)fs{r, —v, T) 
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Here <>cr is used to distinguish the present average of variables that are calculated only at 
the instants of collisions from the preceding ones. Further the integrands in both numerator 
and denominator take into account the effects of particles whose velocities are changed from 
V/ to V at time t, before which a free flight time of r takes place. Thus, fF{r,v/, T,t) is the 
density of particles with velocity v/, before collision, at the point r and time t, which had 
a time of free flight time of length r before the collision, and Ucr is the average collision 
frequency for trajectories on the repeller. Here 

vi = V — 2{v ■ n)n. (25) 

and fi is the unit vector pointing from the center of the sphere (or disk) to the point of 
impact. The factor \h ■ v\ appears in the integrals when one takes into account the rate at 
which collisions take place between the moving particle and the scatterers. Eq. ( P^ should 
be compared with Eq. (^). The mathematical meanings of these averages are, of course, 
different. In Eq. (^ the function to be averaged depends on the random matrix p and the 
phase space variables r and v. In Eq. (|^), the average depends on the random variables 
time of free flight r, collision vector n and phase space variables r and v. But in both cases 
the average is restricted to the repeller. 

The distribution of times of free flight /(r, r, t) and therefore the average on the repeller 
can be expressed as an expansion in gradients of the density. 

/^(f,t;,r,t) = 5{v-Vo){iJo{'r)nm{r,t)+i,({T)v- Vn^{r,t) + ^^{Tyv''n^{r,t)). (26) 

The functions ipf are calculated in appendix A for two and three dimensions. By substituting 
Eqs. ( p6D and (|I^ for fp and fs, respectively, one finds that through second order in the 
gradient, Eq. (E3) assumes the form 



/ dTdn\n ■ v\Q{n ■ v){ij^ - (^^f (1 - 2{n ■ vf) + V^f )g^t;g) In \U ■ e\ 
/ dn\h ■ v\<d{n ■ v){l - -f{l - 2{n ■ vy)q^vl) 

Here a number of simplifications have already been made. All terms that are odd in the 

velocity have been left out since they yield zero on integration. Second order terms in the 
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denominator that give zero on integration over the velocity due to the normahzation of 
the Chapman- Enskog solution have been left out likewise. Furthermore terms of the form 
vv -.VnSJnfiy) have been replaced by ^| Vnp/(f ), which gives the same result on integrating 
over velocity. Finally the factor l — 2v-n results from integrating v' ■ Vn = {v — 2{v- n)n) ■ Vn 
over n. The components of h normal to v give vanishing contributions. By defining the 
constant A = J (in|n-u|(l — 2(n ■ v)"^)/ J (in|n ■ -010(72 ■ v), which vanishes in d = 3 and equals 
-1/3 in d = 2, and expanding Eq. (pTf) up to order g^, we obtain the final result 



Amax(7?.) = i^cr / dTdh\n ■ v\Q{n ■ v) 1^, 



V(^) - + (1 - 2(n ■ vf)) fvl \ (28) 



ln|U 


■ e 






\h ■ 


■ v\ 


|0( 


n ■ v) 



IV. THE LYAPUNOV EXPONENT IN 2 DIMENSIONS 

There is only one positive Lyapunov exponent in two dimensions. Therefore the positive 
Lyapunov exponent can either be calculated with the help of Eq. (plD, i.e. formally as the 
sum of the positive Lyapunov exponents or as the maximum Lyapunov exponent with the 
help of Eq. (|2^). We will use both methods, to demonstrate that the two approaches lead 
to the same result. For the two dimensional Lorentz gas, the radius of curvature matrix 
reduces to a scalar, p, the radius of curvature of two nearby trajectories (see and I), and 
the ELBE has as variables r*, and p. In this case, the ELBE, Eq. (^) has the form ||^J^ 

OF ^ dF d ^ 
-^+v- — +v—F + uF 
at or op 

V /"""/^ a 
= — d(j) cos (j) 6{p COS0) X 

2 J~n/2 2 

fBif,v',t). (29) 

Here we have supposed that the typical value of p before collision is of the order of the 
mean free path length, and the delta function appearing in Eq. (^ has been simplified to 
that appearing in Eq. (^). The solution of the ELBE is normalized according to Eq. (^. 
With if we denote the velocity of a particle immediately before the collision, which results 
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in a velocity v after the scattering event (]25|). For the solution of Eq. (pQ]) we impose the 



boundary condition that F vanishes both at zero and at infinite values of p, i.e., 



F(r, V, p = 0,t) = lim F(r, v, p, t) = 0. 



(30) 



We impose the boundary condition at p = 0, because the radius of curvature increases during 
free motion of the particle, collisions never reduce its value to zero, and all trajectories with 
an initial negative value for the radius of curvature will acquire positive radius of curvature 
with probability one. The negative sign of the radius of curvature is conserved only on the 
stable manifold, which has Lebesgue measure zero in the phase space. 

The integration over (j) can now be performed and we obtain from Eq. (pQ]) 



where a = 2p/a < 1. The velocities v/± result from the evaluation of the delta function 
at precoUisional velocities which satisfy the relation v ■ n = av. We may use the expansion 
Eq. ( p!4D for the Lorentz-Boltzmann density fBirjVyt), by replacing —v in Eq. ([1^) by v. 



Using the Chapman- Enskog Ansatz, Eq. (|T2]), for F, we obtain equations for ipo,ipi and 
1^2 by comparing equal orders of the gradients of the density rim on both sides of Eq. ( p2D . 
The zeroth order equation is ( and I) 




(31) 



ifBir,^f^,t) + fBif,^r_,t)). 



Then we use the relation n{—(j)) = —n{(f)) + 2cos0 v to obtain, from Eq. (|3T|), an equation 
for the distribution function -F(r, v, o", t) = aF{f, v, p, t)/2 given by 




2v d 
a da 



ijQ + uipo = z/e(l - a) 



(1-^2)1/2' 



cr 



(33) 
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with solution 

qi^^-vaa/{2v) for a > 1 

^o(a) = <i (34) 
f [1-(1-(t2)i/2]^ fora< 1 

This solution is continuous at a = 1 and fulfills the normalization condition Eq. up to 
corrections of relative order h = no?. 

By equating the terms to first order in the gradient in the density in Eq. (P^D, we obtain 
an equation for ipi, given by 

+ ?l^''+---!^(i-^)fr^ (35) 

Here we used the fact that the time derivative is of second order in the gradient - via 



the diffusion equation. The solution of Eq. ([35| ) satisfies the normalization condition Eq. 
(p!6|) with C2 = — This can easily be seen by integrating both sides of Eq. (|35| ) with 
respect to cr. The solution of Eq. ( pS]) has to be continuous at cr = 1 and it is given by 

-^(Te-'''^"/(2^.) + iLg-^WCs*;) for (T > 1 
Ma) = { (36) 

(72)1/2(1 + 2^2)], for (7 < 1. 

By comparing the terms in order V^n^ in Eq. (|32|), we obtain the equation for ip2 given 

by 

^M^) + Im^) + ^M^) + —^^2{<y) = (37) 

t>^ 2 a da 

Here we used the diffusion equation for the density rim 

^n^(f,t) = DV2n^(f,t), (38) 

and the Chapman-Enskog solvability condition, Eq. ([17|) . The low density value of the 
diffusion coefficient D in two dimensions is D = (3i;2)/(8z/). For our purposes it is sufficient 
to write the solution of Eq. (|37D in the form 



r<y ' D 1 

M^) = -e~— da'e—i-Mcr') + 7;M'^))- 39 
Jo V 2 
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Once we know ipo, ipi and V'2 we can obtain averages on the repeller using Eq. ([T9|). The 
positive Lyapunov exponent is (see Eq. (|2T|)) 

A+(7^) = '^{-)nep (40) 
a a 

{-)Rev = I M^O + i^^Mo + ^) - i^2yof) - (41) 



Using Eqs. (|3J, ^ |3^) we can now calculate the result in the lowest order in the density 
by standard integrations, and we obtain, to second order in the gradients, 

A+(7^) = A+ + - Df. (42) 

Here \q is the equilibrium solution for the Lyapunov exponent of an infinite system ( @], I) 
given by 

A+ = 2nav[~ \n{2na'^) + 1 - C], (43) 

where C is Euler's constant. 

We can also use the method for the calculation of the largest Lyapunov exponent to 
obtain the result for as an average over the distribution of time of free flight on the 
repeller. The rank of the matrix U, defined in equation (|22|) is unity in two dimensions, i.e 
f/ is a scalar given by 

U = {1 + ^) (44) 
acos(p 

The leading contribution in the low density approximation is obtained by keeping only the 
term proportional to r in Eq. (^4|). The maximum Lyapunov exponent is then in leading 
order 

1 2v 

A+(7^) = -— (ln(— r) - In cos 0),,. (45) 

{r)cr a 



By using Eqs. (|26|) , (pq), and (|A3|) for the distribution of times of free flight in two 
dimensions, we can easily recover the result, Eq. (|42|) . 



17 



V. THE LYAPUNOV EXPONENTS IN 3 DIMENSIONS 



The ELBE for the three dimensional Lorentz gas can be solved by choosing an appropriate 
parameterization of the (ROC) matrix p. Starting from Eq. (1) we can simplify the delta 
function appearing in the restituting part of the collision operator (see I), so as to obtain 
OF ^ dF ( d d \ ^ 

—uF-l — / d(f) da sin (j) COS (f) Y\_ ^ {Pij ~ Pi ^ 

71 JO JO .^^ 

J dp[, J dp[, J dp'^,F{f, v', p', t). (46) 

Here the average (infinite system and low density) collision frequency, is v = na'^vir. 
The ROC matrix is symmetric, therefore only three parameters are necessary to describe its 
behavior. As shown in I it is convenient to use the eigenvalues pi,p2 and the off diagonal 
element pu as parameters. Together with the normalization condition Eq. (P) and the 
solution of the Lorentz-Boltzmann equation as a gradient expansion, given by Eq. (p!4|), 
with —V replaced by v, and C3 = — 1/z^, 

fsif, V, t) =Af6{v- Vo){nUr, t) - ■ Vn„(f, t) + ■■■), (47) 

with A/" defined below Eq. We find that Eq. (^) can be written as 

— + v — + v(— + —\f + uF- 
dt df \ dpi dp2 J 



— I d(j) I (iasin0cos0 I cos2a| 5(pi cos0)(5(p2 

TT JO JO 2 



2 COS0 

^(Pi2 + (acos0 tan^ sin 2a;) /4) fB{r,v',t). (48) 
The factor |cos2a| results from the transformation from the variables pii,P22,Pi2 to 



Pi^ P2j Pi2 in the 5-functions in the right hand side of Eq. (^), and we set -F(pi, P2, P12) = 
F'ipn, P22, Pu)- Please note, that F is not yet normalized with respect to pi,p2,Pi2- The 
velocity v/ before the collision depends on the collision angles a and (p, through the relation 

vf = V — 2{h ■ v)h, with 

n = cos -y + sin cos a Vj_^i + sin sin a Vj_^2, (49) 



and the unit vectors v,v±^i,v±^2 form an ortho-normal set. 

Now, using Eq. (0), we can perform the a integration, and we keep only the non- 
vanishing terms in the integration. After integrating over a and making fs explicit, we 
find that we can rewrite Eq . (|48|) as 



— +v —+v (— + —\ F + uF 
dt dr \ dpi dp2 1 



vra Jo 



sin d) cot (hQ [l 



4p 
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acos(/)tan^ - 



a cos I 



2 cos I 



-)5{v - vq) (nm{r,t) - -(1 - 2cos^0)tT- Vnm{r,t) + 



(50) 



Here we have explicitly indicated only the terms needed for our further calculations. It is 
now convenient to introduce the definitions 2pi/a = ai and 2pi2/a = ai2, and to perform 
the (J) integration. We find that 

dF ^ OF , f d d \ ~ ~ 



e 1 



2criO-i2 



\-ai 



1 - (Tf (Ti 



5(f - Vq) {nm{r,t) - -^(1 - 2o\)v- Vnm{r,t) H ^ 



(51) 



We point out that F is not normalized with respect to 0"i,ct2 and ai2, since we have 
not yet introduced the appropriate Jacobian. Equation ([51| ) can be further simplified by 
introducing a new set of variables (0<s,0<z<l, — 7r/2 < 7 < 7r/2), defined by the 
relations 



CTi = z + s. 



0-2 = - + S 

z 

(Ji2 = sin(7) 



(52) 



2z 



The normalized distribution function f{f,v,s,z,'j,t) 
obeys the equation 

-^ + df + ——f + vf 
ot a OS 



J (s, z, -f)F{r, V, pi, p2, P12, t) 
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'^^^-z 6(1 - z)Q{z)Q{t:/2 - 7)6(7 + 7r/2)(5(s - 0+) (53) 



5{v- Vq) [nm{r, t) - -^(1 - 2z^)(9n„(f, t)^ , 



where d = v - V , and J(s, 7) is the Jacobian given by 

j( X |5(Pll,P22,Pl2)|5(Pl,P2,Pl2) ||5(C^1,C^2,C^12)| /^.x 

J(2;,s,7)= — -|— \\\—^( ^ V>V 

o{pi,p2,pi2) o{ai,a2,ai2) d{s,z,-f) 



a 



3 ^l + z^){l-z^ 



16 z 



3 



(55) 



We point out that this Jacobian is independent of s. It guarantees the proper normahzation 
of / as 



roc rl _ 

/ ds dz d-ff{f,v,s,z,-f,t) = fB{r,v,t). (56) 

Jo Jo J~n/2 



where fs is the solution of the standard nonequihbrium Lorentz-Boltzmann equation. 
The physical meanings of s, z and 7 become a bit more transparent if we consider Eq. 

(0) and Eq. (^) in more detail. We begin by noting, that the distribution function, 
F{f, V, pi, P2, P12, t), for the two eigenvalues and the off diagonal element of the ROC matrix 
is established through the dynamic process involving intervals of free flight separated by 
collisions of the moving particles with the scatterers. The time dependence of the ROC 
matrix can be completely expressed in terms of scattering angles and time of free flight. In 

(1) it was shown that the eigenvalues of the ROC matrix increase linearly in time during 
a free flight and the off diagonal element stays constant. The reparameterization Eq. (|52D 
exactly reflects this behavior. We can see that the dimensionless parameter s corresponds to 
the time of flight, which is also clear since the 6 - function in s in the right hand side of Eq. 
(^) shows that the gain term is a source of particles with a free flight time of zero. From 
Eq. (|52|) we can then conclude that z, respectively ^, correspond to the possible values of 
normalized eigenvalues cti, cr2 at the collision. Equation (|48|) shows that z is physically the 
cosine of the scattering angle (p. With the same arguments, it can be seen that up to a factor 
of 2, 7 finds it physical correspondence in the azimuthal scattering angle a. 
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As one might suspect here, the solution of Eq. ( p3D can be interpreted as a joint distri- 
bution function for time of free flight and collision parameter. This may seem surprising at 
first, because the latter are statistically independent quantities. Correlations in / however, 
are the result of considering this function at fixed time and position. In a comoving frame, 
i.e. cconsidering /(r — ^av,v, s, Zyjyt — aa/2v), one would find the variables s, z and 7 to 
be uncorrelated indeed. 

The distribution of times of free flight in 3 dimensions Eq. (^l|) can be recovered, if we 
integrate Eq. (|53D over z and 7 and identify as/{2v) with r. 

(2 rT^/'i rl _ 2,V 

fF{f,v,T) = — d-f dzf{f,v,—T,z,-f) (57) 

2v J-n/2 Jo a 

For large systems, we can solve Eq. (0), as in the 2 dimensional case, by using a gradient 
expansion 

/ = 6{v - Vo)Afe{n/2 - 7)0(7 + n/2){n^{f, t) z) + ^dn^{r, t) ^i(s, z) 

2v 

a? 

+—V^nra{f,t)ilj2{s,z) + ■■■). (58) 

The solution of Eq. (|53|) has to satisfy the relation (|56|) The quantity ?/'o('S, z) was already 
obtained in I. It satisfies the equation 

-^V-o + z>^o = — 0(1 - z)Q{z)z5{s - 0+), (59) 
as 71 

where the dimensionless collision frequency is u = ^u. The solution of Eq. ( ^9]) is 

^^0(5, z) = — 0(1 - z)Q{z)z Q{s)e-^'. (60) 

TT 

The equation for ipi is obtained by keeping only the terms proportional to drim 

i'o + ^^1 + z>V^i = — 0(1 - 2)0(2)2(1 - 22^) 6{s - 0+). (61) 
as 7r 

This equation is easily solved to give 

^i(s, z) = ^ {sz + ^2(1 - 22^)) 0(1 - z) 0(2) 0(s)e-^^ (62) 

TT V 
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Considering the terms of order V^nm(^, t) and keeping only the scalar part of dd, we obtain 
the equation for 'ip2 

Di^o + li'i + + = 0. (63) 
3 as 



Here the Chapman Enskog solvability condition Eq. (|T^ and the diffusion equation for n^, 
Eq. (|38|) were used. We also introduced the dimensionless diffusion coefficient D = —D. 
The solution of this equation is given by 

^2(s, z) = — {-bsz + + ^sz{l - 2^2)) 0(1 _ z) Q{z) Q{s)e-^\ (64) 

TT b OP 

The formulas for the sum of the Lyapunov exponents and the maximum Lyapunov expo- 
nent in equilibrium were obtained in I. The same formulae are valid in the nonequilibrium 
case, if the averages are replaced by averages on the repeller. For the sum of the Lyapunov 
exponents Eq. (|19D gives 

f A,(7^) = ((^ + —^))Rer>, (65) 
2v z + s 1/z + s 

(66) 

while Eq. (^) leads to 

XmaxCR-) = z/cr(ln|U(cos0, a,r) -e^Dcr, (67) 
with the average collision frequency on the repeller Ucr = -j^^ and 



\V{z,a,T)-e^\ = —]J — (68) 

with z = COS0, as derived in I. The angle ip specifies the direction of the unit vector in the 
plane perpendicular to the trajectory. In general an additional average over its stationary 
distribution P{tp) is necessary. However, it can easily be shown that the corrections to an 
isotropic distribution of directions P{ip) are at most of order cos('?/')V^n(r, t) and do not 
contribute to the average in Eq. (^) in order V^n(r, t) (see appendix B). 
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Now, using the equations Eqs. |T9D, we are led to the determination of the sum of 
the positive Lyapunov exponents, and we obtain, to second order in the gradients, 

A+„.(7^) + Kmnm = his - Df (1 + 2(ln(n/2) + C)), (69) 

with h = nna^ = 2u. Here, h'j^g is the KS-entropy for an infinite system [I] at equihbrium 
given by (|70|), C is Euler's constant, and q is defined by (|20|), and 

/i?^5 = 2na\7r[-ln(n/2)-C] + ---. (70) 

with the dots indicating higher density corrections. 

This expression for the sum of the Lyapunov exponents can also be calculated by aver- 
aging In det U = Trace In U over the distribution of times of free flight and the scattering 
angles Q with the help of Eq. ( ^8]) , and by replacing In | C7 ■ e| with Trace In U In leading 
order for large times of free flight i.e. small density of scatterers, we obtain 

+ A+.J7^) = z/,,(Traceln(^p;V)),,. (71) 

The matrix p_^_ is the ROC immediately after a scattering event. It depends in leading order 
in the density only on the scattering angles (p and a and is defined in paper I. Here it is only 
important to notice that it has eigenvalues pi = |cos0 and p2 = ^ ^ . The trace of the 



logarithm in Eq. (71) leads therefore to a cancellation of the terms depending on In cos 0, 



so that we only have to evaluate 

2v 

K.a.m + A+.J7^) = 2z/,,(ln(-r)),,. (72) 

a 

This strategy also leads to the result Eq. ([691). 

With the help of Eqs. (^ ^ |2^) and Appendix A the maximum Lyapunov exponent 



is given, to second order in the gradients, by 

AL.(^) = Al, + Dq\- ln(n/2) - C + ^ - In 2) , (73) 

where X^^ax is the equilibrium value of the maximum Lyapunov exponent for an infinite 

system 
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Al. = na'vT^i- ln(n/2) + In 2 - 1 - C] + ■ ■ • . (74) 

The expression for the smallest positive Lyapunov exponent can be obtained, to second 
order, from Eqs. (|69D and ([73|), as 



A^.„(7^) = AL„ + Dt{- ln(n/2) - C - ^ + ln2) ■ ■ • , (75) 

with 

ALn = ^«'^^[- ln(n/2) - ln2 + 1 - C] + ■ ■ • . (76) 



VI. DISCUSSION 

We have now calculated, analytically, the spectrum of positive Lyapunov exponents on 
the repeller for the open, dilute, random Lorentz gas in two and three dimensions. Then 
using the escape-rate formula we may infer the values of the KS entropies on the repeller 
as well. We find that the corrections to the equilibrium values of the Lyapunov exponents 
and KS entropies are of order where L is some characteristic size of the open system. 
We should point out that Gaspard has discussed a reformulation of the escape-rate formula 
so as to be able to express the diffusion coefficient in terms of the Hausdorff dimension of 
the fractal repeller 0. Using his method we can easily see that the dimension of the fractal 
repeller is slightly less than the embedding dimension (3 for d = 2, and 5 for d = 3) by 
terms of order These results are to be expected for the Lorentz gas, since diffusion is 

normal, and the escape-rate formula should be free of difficulties. 



Gaspard and Baras [|rT| have examined the chaotic properties of the periodic open Lorentz 
gas at sufficiently high densities that there are no infinite horizons for the moving particle. 
They used numerical simulations, and obtained independent results for the Lyapunov ex- 
ponents and for the KS entropy on the repeller, as functions of the system size, L. Then 
they compared these results with numerical and with approximate analytical values for the 
diffusion coefficient and found good agreement. 
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At the present time there are no computer simulations of open, random Lorentz gases, 
to which our analytic results can be compared, but we may compare our results to intuitive 
expectations. To do this we will need to note that the average collision frequency, Ucr or 
equivalently, the mean free time between collisions on the repeller, = l/t'cr, differ from 
the corresponding quantities in an infinite system. In fact r^r can explicitely be calculated 
by using an expression analogous to Eq. (^) , where In | U ■ e | is replaced by r and the 
factor Ucr is dropped. Together with the expressions for ipf in two and three dimension Eqs. 
( p6| , |A^ , and we obtain an expression for the mean free time between collisions on the 
repeller given by 

Tcr = — = ^g'. 77 

This result shows that the mean free time on the repeller is smaller than that in the 
infinite system. In addition, its dependence on the escape-rate and equilibrium collision 
frequency is the same in two and three dimensions. Thus, we see that trajectories on 
the repeller are constrained to have higher collision frequency than those for the infinite, 
equilibrium system. 

The Lyapunov exponents i.e. the rates of separation of trajectories, are larger on the 
repeller than in the infinite system. One might have expected, that due to the restriction of 
the available phase space volume for a particle on the repeller, the rate of separation might 
have been smaller on the repeller, but the increased scattering rate counteracts this effect. 
However, if we compensate for this effect by expressing the Lyapunov exponents in units of 
the mean free time on the repeller, the terms proportional to lead to a decrease of the 
Lyapunov exponents because on average ln(A/a) decreases. Using Eqs. ( ^21) and (^) for two 
dimensions and Eqs. ( |69| , [73| , |75D and ( |77D for three dimensions, the Lyapunov exponents in 
natural units on the repeller are given by 

^+ Y Y) 

\'^{JV) Tcr = — — - — q"^ in 2 dimensions, (78) 
and for three dimensions 
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V u \ 4 J 
A+ (n)T -A° ^-2/ 7-81n2 \ 

(AL.(^) + >^Lnm) = % - (79) 

It is also interesting to compare the KS entropy on the repeller with its value in an infinite 
system. To make this a bit more concrete, we consider a slab geometry, i.e. absorbing walls 

at X = ±L/2, and 

Q = l- (80) 

Then with Eqs. ( ^21) and (p9| ) respectively, we obtain. 

h° 3 

/^^s(7^) = + ~ 2'* ^ (^/^)^ two dimensions (81) 
/^^g(7^) = h^^s + - 2) D{'n/Lf in three dimensions (82) 

Thus the KS entropy increases above its infinite system value when measured in standard 
time units. As in the case of the Lyapunov exponent, this trend is only due to the increased 
scattering rate on the repeller. When measured in natural units on the repeller, the KS 
entropy has an especially simple form. 

3D 

hKsiT^) Tcr = — ^ — i'^/^y ill t'^o dimensions (83) 

u 2 u 

h° D 

hKsiT^) Tcr = - 2 —(n/Ly in three dimensions (84) 
u V 

We conclude with a number of points: 

1. The principal problem with the escape-rate method as an analytical method for com- 
puting transport coefficients, apart from the inherent difficulties involved in calculating 
dynamical quantities, is that as yet we have no analytical methods for calculating the 
KS entropy on the repeller, independently of the escape-rate formula. At the mo- 
ment we can only use analytic techniques to calculate the diffusion coefficients and the 
Lyapunov exponents, leaving the KS entropy as a quantity to be derived from them. 
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2. It would be valuable to have some results from computer simulations with which to 
compare the results obtained here. 

3. As mentioned earlier, the information and Hausdorff dimension of the fractal repellers 
in both two and three dimensions are very close to the full phase space dimensions, 
three, for the Lorentz gas on the plane, and five, for the Lorentz gas in space, but are 
smaller than these values by terms of order . This is a simple consequence of the 
application of Kaplan- Yorke type formulae to fractal repellers [0 . For the information 
dimension we obtain 

= 3 - 2-^ + 0(1/L^) forrf = 2 (85) 
rf^ = 5 - 2-^ + 0(1/L^) forrf = 3 (86) 

where 7 is the escape rate and Ag , A^j„ are the positive Lyapunov exponent, and the 
smaller of the positive Lyapunov exponents in the infinite system, respectively. Gas- 
pard and Baras have used these fractal dimensions to express the diffusion coefficient 
in terms of the Hausdorff and information dimensions of the fractal repeller for the 
two dimensional case |]6|,|TT|. 



4. A problem for further study is to extend the calculations given here to a system of 
many interacting particles such as gases of hard disks or hard spheres. Some progress 
in this direction has been made, and it is now possible to get analytic results for the KS 
entropy and the largest Lyapunov exponents for dilute hard disk or hard sphere gases, 
in equilibrium, in the thermodynamic limit |]T^,|r^,|T^ . It would be very interesting 
to apply the escape-rate formalism to the transport coefficients, such as the shear 
and bulk viscosities, and thermal conductivity, appropriate for fluid systems, using 
the method of Ref. P], and to determine the effects of the fractal repeller on the 
dynamical quantities. 

5. The thermodynamic formalism for hyperbolic chaotic systems provides a very useful 
method for expressing many of the chaotic properties of both open and closed systems 
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in terms of one quantity, the topological pressure 0. The use of kinetic theory to 
evaluate the topological pressure for a dilute random Lorentz gas should certainly be 
possible, but has not yet been undertaken. 

6. In the next paper in this series we will consider the case of a Lorentz gas with a 
charged moving particle placed in an external electric field as well as in a random 
array of scatterers. Then a Gaussian thermostat is applied which keeps the kinetic 
energy of the moving particle fixed. The system eventually reaches a nonequilibrium 
steady state. We will calculate dynamical properties of the moving particle in this 
nonequilibrium steady state and compare with the results of computer simulations. A 
preliminary version of this work has already been published |19,20t|. 
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APPENDIX A: THE DISTRIBUTION OF TIMES OF FREE FLIGHT 

It is useful to have an expression for the mean collision frequency or, equivalently,for the 
mean free time between collisions for trajectories on the repeller. Here we argue that such 
expressions can be obtained in a simple way from a Lorentz-Boltzmann equation, modified so 
as to include a new variable r which is the time since the last collision. For the distribution 
of particles surviving in the system at time t, with position r, velocity v, and with time r 
since the last collision we propose the equation 

d ■ d d 

—Jf + V ■ Vfp + V ■ —fp + -rrfp = -vfp 

Ot OV OT 
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(Al) 

with 

fp{f,v',T = 0,t) = na'^-^v J dne{-n-v)\h-v\fB{f,v',t). (A2) 

Eq. (^T|) can be derived in a way similar to the heuristic derivation of the usual Lorentz- 
Boltzmann equation, but a few changes have to be made. Scattering and absorption at the 
boundary are the only mechanisms which will reduce the number of particles with a certain 
time of free flight. Scattering will also act as a source of particles but always with a time 
of free flight r = 0. This effect is taken care of by the initial condition Eq. ( |A2|) . The 
dependence of the distribution of times of free flight on both position and velocity is due to 
the non-uniformity in these variables of Jb resulting from the absorbing boundary condition 
in combination with a higher survival rate for particles that collide more frequently (they 
diffuse more slowly). We note one important difference, which will become crucial, if the 
effect on external fields is considered. In the Boltzmann equation for the one particle phase 
space density J'b, which is based on the consideration of the number of particles in a fixed 
volume element in phase space, the streaming term is derived from t)) = with 

q = (r, v). This is due to the conservation of particles in the absence of scattering events. In 
the derivation of fp we have to consider the number of particles with a certain times of free 
flight T after the last scattering event and count how many of them are still there a time step 
dt later. We therfore have to use a comoving frame, which means that the streaming part 
of the equation for fp has the form q-^friQ, t). A term analogous to cannot appear in 
a comoving frame since it counts the difference of ingoing and outgoing particles in a fixed 
phase space volume element. 

The solution of this equation can be obtained in two and three dimensions as a gradient 
expansion (see Eq. ^6]). The solution strategy is completely analogous to that used for 
solving the ELBE in two and three dimensions. For two dimensions we obtain 

V^o^(r) = z/e— 
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V^f (r) = (- - vr)e-'^^ 

VAr) = \{-\yr ^\{yrf)e-^\ (A3) 
In 3 dimensions we obtain 

(r) = -z/re""" 

^r(r) = i-(-z.r + i(z.r)2)e— . (A4) 

As stated before above Eq. ^ the same results may be obtained by integrating / over 
the variables z and 7. 



APPENDIX B: THE DISTRIBUTION OF EIGENDIRECTIONS P(V') 



In order to calculate the maximum Lyapunov exponent of a product of uncorrelated 
random 2x2 matrices n^i(T, 0, «), we must take into account the fact that the matrices 
Uj do not in general commute with each other. In order to use standard theorems to 
calculate the largest eigenvalue of a product of random matrices, we have to determine the 
distribution P(V^) of the angle V', generated by acting with the random matrix U on the 
unit vector e(-?/') = (cos('?/^), sin('?/')) |T^. If this distribution is not isotropic in t/^, the proper 
form of P(V') must be determined from the solution of an appropriate Frobenius-Perron 
equation. 

In our case, the Frobenius-Perron equation for this distribution is given by 



2-K 



Jo 



(Bl) 



The angle 4>, a) is implicitly defined by 

U • e^' = |U ■ e^/ 



^ / \ 

cos tpi 

sin-j/'i 



(B2) 
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where 0, a are the scattering angles. The distribution function has to obey the normahzation 
condition 



/•27r 

/ d^PH)) = 1 (B3) 
Jo 

The average in Eq. ( [Bl| ) (...) is in our case the average over the distribution of times of 
free flight, scattering angles, velocities and space multiplied with the survival probability, 
as defined in Eq. (^). Since, contrary to the case of motion in an external field ||20|, the 
equation of motion for U is the same for an open system as for an infinite system when there 
is no external field, we can use here the expression of |U ■ e^| derived in I. To evaluate the 
ensemble average in the Eq. (|B1|), we have to use the gradient expansions for the survival 
probability and the distribution of times of free flight Eqs. ([T^ , p!2D . Since the ensemble 
average involves an average over the velocities, the term proportional to vVn and all other 
non scalar terms vanish. Therefore the deviation of P{ip) from its isotropic value, {2tt)^^ can 
at most be of order V^ra. Then, due to the normalization condition Eq. ( P3| ) this deviation 
must also be be proportional to cos{milj), where m is a positive integer. Consequently the 
only additional term of order V^n in the expression for the maximum Lyapunov exponent 
would come from multiplying the zeroth order term (in the gradient of n) of ((In |U ■ e('?/^)|) 
(which is independent of ip) with the term of order V^n in the expression for P{ip). But 
this averages to zero due to the factor cos(m-?/'). Thus we may use an isotropic distribution 
in the angle if) when calculating the largest eigenvalue of the product of random matrices, 
and we then obtain Eq. (|73D. 

We can understand this result also in a more intuitive way. This demonstration proceeds 
in four steps. 

1) An isotropic precollisional distribution of ip will give rise after scattering to anisotropies 
of order —(f. To understand this one should note that angles if) and ip^n can be identified. 
This merely amounts to an interchange of reference and tangent trajectory. If the contri- 
butions of such angle pairs are averaged, the term linear in v ■ Vn cancels. The sole source 
for the anisotropy in this case is the anisotropy of the survival probability as a function of 
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velocity. 

2) A completely anisotropic precollisional distribution, of form 6{tp — ipo) say, will give rise 
to a postcollisional distribution that is isotropic in tp, up to corrections of order — g^. This 
is a consequence of the isotropic distribution of the azimuthal scattering angle a resulting 
from the random distribution of scatterers. Again, the anisotropy is due to anisotropy of 
the survival probability and the linear term in v ■ Vn vanishes due to identification of ip and 

Ip ±TT. 

3) If we require that the anisotropic part of the distribution of ip cancels on averaging, so 
does the isotropic part of the distribution just after a collision, that is, in an ensemble av- 
erage in which we put together all contributions of ensemble members with a collision at an 
equal position and time. 

4) Now since the precollisional anisotropic distribution itself is of order according to 1) 
(and nothing will be added to that just because of the result emerging here), its contribution 
to the postcollisional anisotropy has to be of order g^. 

The postcollisional anisotropics resulting from the isotropic part of the precollisional distri- 
bution are fully accounted for by Eq. ^ and so this expression gives correct results through 
order g^. 



32 



REFERENCES 

[1] H. van Beijeren, A. Latz, and J. R. Dorfman, Phys. Rev. E, 57, 4077, (1998). 

[2] P. Gaspard and G. Nicolis, Phys. Rev. Lett., 65, 1693, (1990). 

[3] T. Bohr and D. Rand, Physica D, 25, 387, (1987). 

[4] P. Gaspard and J. R. Dorfman, Phys. Rev. E, 52, 3525, (1995). 

[5] J. R. Dorfman and P. Gaspard, Phys. Rev. E, 51, 28, (1995). 

[6] P. Gaspard, Chaos, Scattering, and Statistical Mechanics, Cambridge Univ. Press, Cam- 
bridge, (1998), Ch. 6. 

[7] J. R. Dorfman, An Introduction to Chaos in Nonequilibrium Statistical Mechanics, Cam- 
bridge Univ. Press, Cambridge, (1999). 

[8] H. van Beijeren and J. R. Dorfman, Phys. Rev. Lett., 74, 4412, (1995); Erratum: Phys. 
Rev. Lett., 76, 3238, (1996). 

[9] P. Gaspard and F. Baras, in Microscopic Simulations of Complex Hydrodynamic Phenom- 
ena, M. Mareschal and B. Hohan, eds., NATO AST Series, No. 292, Plenum, New York, 
(1992). 

[10] E. Ott, Chaos in Dynamical Systems, Cambridge Univ. Press, Cambridge, (1992). 
[11] P. Gaspard and F. Baras, Phys. Rev. E, 51, 5332, (1995). 
[12] B. Hunt, E. Ott, and J. Yorke, Phys. Rev. E, 54, 4819, (1996). 

[13] S. Chapman and T. G. Cowhng, Mathematical Theory of Non- Uniform Gases, Cambridge 
Univ. Press, Cambridge, (1970). 

[14] Ya. G. Sinai, Russ. Math. Surv.,25, 137, (1970). 

[15] A. Crisanti, G. Paladin, and A. Vulpiani, Products of Random Matrices in Statistical 
Physics, ( Springer- Verlag, Berlin, 1993). 

33 



[16] R. van Zon, H. van Beijeren, and Ch. Dellago, Phys. Rev. Lett., 80, 2035, (1998). 

[17] H. van Beijeren, J. R. Dorfman, H. A. Posch, and Ch. Dellago, Phys. Rev. E, 56, 5272, 
(1997). 

[18] J. R. Dorfman, A. Latz, H. van Beijeren, CHAOS, 8, 444, (1998). 

[19] H. van Beijeren, J. R. Dorfman, E. G. D. Cohen, H. A. Posch, and Ch. Dellago, Phys. 
Rev. Lett., 77, 1974, (1996). 

[20] A. Latz, H. van Beijeren, and J. R. Dorfman, Phys. Rev. Lett. 78, 207, (1997). 



34 



